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Abstract 

In terms of its highly pathogenic nature, there remains a significant need to further define the immune pathology of SARS- 
coronavirus (SARS-CoV) infection, as well as identify correlates of immunity to help develop vaccines for severe coronaviral 
infections. Here we use a SARS-CoV infection-reinfection ferret model and a functional genomics approach to gain insight 
into SARS immunopathogenesis and to identify correlates of immune protection during SARS-CoV-challenge in ferrets 
previously infected with SARS-CoV or immunized with a SARS virus vaccine. We identified gene expression signatures in the 
lungs of ferrets associated with primary immune responses to SARS-CoV infection and in ferrets that received an identical 
second inoculum. Acute SARS-CoV infection prompted coordinated innate immune responses that were dominated by 
antiviral IFN response gene (IRG) expression. Reinfected ferrets, however, lacked the integrated expression of IRGs that was 
prevalent during acute infection. The expression of specific IRGs was also absent upon challenge in ferrets immunized with 
an inactivated, AI(OH) 3 -adjuvanted whole virus SARS vaccine candidate that protected them against SARS-CoV infection in 
the lungs. Lack of IFN-mediated immune enhancement in infected ferrets that were previously inoculated with, or 
vaccinated against, SARS-CoV revealed 9 IRG correlates of protective immunity. This data provides insight into the 
molecular pathogenesis of SARS-CoV and SARS-like-CoV infections and is an important resource for the development of CoV 
antiviral therapeutics and vaccines. 
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Introduction 

Severe Acute Respiratory Syndrome (SARS) disease hit the 
world in late 2002 and in 4 months swiftly spread to 29 countries 
infecting over 8,000 people and killing over 700 [1]. The 
etiological agent of SARS disease was determined to be of the 
coronavirus (CoV) family; the largest family of single-stranded, 
positive-sense RNA genomes known [1], The overall mortality 
rate of SARS corona virus (SARS-CoV) infection was —10% but 
this rate was 50% in patients over 65. Prior to the emergence of 
the SARS virus, coronaviruses were known to cause mild upper- 
respiratory tract diseases in humans. In contrast, SARS-CoV 
infection caused severe disease in the lower respiratory tract 


disease with symptoms ranging from flu-like and viral pneumonia 
to acute respiratory distress syndrome (ARDS) and fatal outcome 
[2-5]. The virus emerged from the Guangdong Province in China 
where it crossed to humans from a zoonotic reservoir. The most 
established theory puts horseshoe bats as the ultimate reservoir for 
the SARS-CoV and implicates palm civets as the intermediate 
species that passed the virus to humans [1]. Aggressive public 
health intervention strategies are credited with successfully 
minimizing the SARS-CoV infection range, although it is 
uncertain if these same public health strategies would sufficiently 
contain a future SARS-CoV or SARS-like-CoV outbreak due to 
virus evolution. 
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Importantly, coronaviruses have a propensity toward frequent 
host-shifting events and over the past 30 years there have been 
many CoV cross-species transmission incidents giving rise to new 
animal and human CoV -based diseases. Coronaviruses infect a 
broad range of species lending further chance for recombination 
events and the advent of new CoV species. Moreover, 
coronaviruses can change cell type, tissue and host species 
barriers with ease [6,7]. Typically, the spike (S) protein of 
coronaviruses determines the host infectivity and the organiza¬ 
tion of the SARS-CoV S protein shows significant similarity with 
other aggressive class I viral fusion proteins: influenza virus HA, 
HIV-1 Env, Simian virus 5, and Ebola virus Gp2 [1], The 
promiscuity of coronaviruses coupled with the tendency for 
mutations to occur gives reason for concern that another CoV 
outbreak is likely and highlights the need for continuous viral 
surveillance and forward development of CoV vaccination 
strategies and therapeutics. 

Although entry of SARS-CoV into mammalian cells has been 
determined to be facilitated by the angiotensin-1 converting 
enzyme 2 (ACE2) molecule [8], the mechanisms by which the 
virus evades host immune responses causing generalized inflam¬ 
mation, increasing viral burden, and severe lung pathology still 
remain a significant scientific problem. Previous studies have 
shown substantial problems with potential CoV vaccines where 
the vaccines cause disease exacerbation opposed to initiating 
immunological protection [9,10]. Recently, several groups have 
described the immunologic response during SARS-CoV infection 
[11] and some have investigated the use of a mouse adapted 
SARS-CoV in the mouse model [12-15], The mouse-adapted 
SARS-CoV (MA15) is a valuable animal model for investigating 
the immune response and possible therapeutic and prophylactic 
strategies for SARS-CoV disease. Although the model helped to 
elucidate immune-pathological events during SARS-CoV infec¬ 
tion and protection [12-15], the caveat of this model is that it is 
based on an adapted virus and not a wild-type SARS-CoV that 
has naturally occurred in nature and cause disease in humans and 
animals. Although death is not observed in our wt TOR2 SARS- 
CoV ferret model, there are still several advantages where the use 
of both models is perhaps of equal importance as results from the 
mouse model compliment findings from the ferret model and vice 
versa. Specifically, use of the ferret model provides several benefits. 

As mentioned above, ferrets are susceptible to wild-type SARS- 
CoV infection from strains isolated from humans [16,17]. 
Furthermore, when infected with respiratory viruses including 
the SARS-CoV ferrets display many of the symptoms and 
pathological features as seen in infected humans as ferrets and 
humans have similar lung physiology [18-21]. Quantitative 
clinical signs displayed by ferrets include a rise in core body 
temperature (fever), nasal discharge (sneezing and runny nose) and 
weight loss [16], 

Here we investigated the immune response transcriptome of 
SARS-CoV pathogenesis in a ferret model infected with an 
unadapted SARS-CoV and subsequently evaluated gene expres¬ 
sion signatures induced with SARS-CoV reinfection. Further¬ 
more, ferrets were immunized with a SARS-CoV vaccine and 
then challenged to compare immunological profiles with the 
SARS-CoV reinfected animals. The objective of this study was to 
identify immune correlates of protection upon reinfection with 
SARS-CoV in ferrets and provide a comprehensive profile of an 
effective and nonpathological immune response to SARS-CoV 
challenge following immunization. This information will not only 
provide a foundation for direct comparison with future SARS 
vaccine studies, but will also allow us to determine what immune 
mediators are responsible for the successful antiviral response. 
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Results 

Effective Immune Responses to SARS-CoV Reinfection 

Previously, the lack of a representative SARS-CoV infection 
animal model has limited the ability to uncover immunopatho- 
genic mechanisms of SARS and has impeded progress in 
vaccination strategies. Currently, a mouse model for SARS-CoV 
infection has described aspects of SARS-CoV pathogenesis 
although these studies utilized a mouse adapted SARS virus 
[12,13,15]. The ferret, Mustela putorius jiiro, displays many of the 
symptoms and pathological features seen in SARS-CoV infected 
humans and is susceptible to unadatped SARS-CoV strains 
therefore suggesting it as a useful animal model for the study of 
SARS-CoV infection and vaccination strategies [5,22—24], Here 
ferrets were assigned to three treatment groups for this study: (A) a 
mock infection control group, (B) a single infection group, and (C) 
a group that was infected at the same time as group B, then 
subsequently reinfected 29 days later. The infected groups were 
innoculated intranasally (IN) with 10' 5 TCID 50 (50% tissue culture 
infective dose) of SARS-CoV TOR2 strain, while the mock 
infected animals received intranasal instillation of vehicle (serum- 
free cell culture medium (SFM)). 

Ferrets were monitored for clinical signs of disease twice daily. 
The most significant symptom observed in SARS-CoV infected 
animals was sneezing, which was observed in one third of the 
infected ferrets between 4 and 11 days post-infection (DPI). 
Sneezing was not observed in ferrets from the control group (A), 
nor in group (C) at any time post- reinfection. 

Neutralising antibody titers rose above baseline from 6 DPI to 
maximal post-infection levels at 15 DPI (Figure 1). Specifically, 
neutralizing antibodies in the single infection group (B) increased 
from study Day 8 to a peak titer (820U) on Day 15. Neutralizing 
antibody titers fell slightly to 720U on Day 30 (day of 
reinfection). Increases in neutralizing antibodies directly corre¬ 
lated with decreased viral burden in the lung (Figure 2). At four 
weeks post-infection, animals in group C were reinfected with 
SARS-CoV and a second spike in neutralizing antibodies 
(1300U) above levels seen during the infection phase was 
observed within a week. Neutralizing antibody titers remained 
elevated for approximately two weeks, then diminished to 400U 
by the end of the study (day 58). The boost in antibody titer may 
have contributed to the lower viral burden observed during the 
reinfection phase (Days 32 to 58). This second peak in 
neutralizing antibody titres was coincident with restricted viral 
replication (Figure 2). Neutralizing antibodies were not detected 
in the mock-infection control group. The TOR2 strain of SARS- 
CoV replicated well in the upper respiratory tract (data not 
shown) and lung within the first week post-infection. Significant 
viral titers above background were detected in lung, ranging 
from one to three log above background. Peak viral titers were 
observed between 3 to 6 DPI in the lung (Figure 2). Viral titers 
post reinfection (Days 30 to 58) were generally restricted below 
2 log values (Figure 2). Furthermore we investigated the lung 
histopathology following SARS-CoV infection and following 
SARS-CoV infection-reinfection (Figure 3). Ferret lung sections 
were obtained front multiple lung lobes at 7 d post primary 
infection or 7 d post secondary infection and stained by 
hematoxylin and eosin. Primary SARS-CoV infection produced 
inflammation and the appearance of lung immune cells primarily 
surrounding small-to-medium bronchial airways at 7 days 
following challenge (Figure 3 B). Ferrets that received infec¬ 
tion-reinfection were largely protected front lung histopathology 
(Figure 3 C). Taken together the ferret mounted an effective 
immune response to the reinfection with SARS-CoV, with 
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Figure 1. Mean levels of serum neutralizing antibody to SARS-CoV. Neutralizing antibody levels in serum were determined as outlined in the 
materials and methods section. Inverse neutralization titre is represented on the y-axis vs. day of study on the x-axis. Note that ferrets were infected 
on study day 1 and reinfected ferrets were also innoculated on study day 30. Mock infected animals received an intranasal instillation of serum-free 
media on study day 1 and had undetectable titres. Values shown represent group mean of 3-4 ferrets per group, and error bars show standard 
deviation. 

doi:10.1371/journal.pone.0045842.g001 


increased neutralizing antibody titre, restricted virus replication, 
diminished clinical symptoms and less lung pathology in 
comparison to the initial infection of the naive host. 

Lack of IRG Expression in SARS-CoV-reinfected Ferrets 

To characterize host immune responses in an infection/ 
reinfection ferret model of SARS-CoV pathogenesis, we 
performed gene expression profiling on lung tissue taken at 
necropsy at days 2, 3, 5, 7, 14 and 28 days post-infection (DPI) 
and days 2, 5, 7, 14 and 28 days post-reinfection (DPR). Three 
ferrets at each time point were profiled and gene expression 
signatures associated with pulmonary immunopathology during 
SARS were identified using an Extraction and Analysis of 
Differential Gene Expression (EDGE) differential time course 
microarray analysis [25] on the SARS-CoV infected ferret lung 
samples. A heat map overview containing 3454 genes identified 
as significantly differentially expressed across all time points and 
classified by Ingenuity Pathway Analysis (IPA) into three clusters 
of functionally related genes (Figure 4); an IL-6/complement, 
an IRG, and an adaptive immune gene cluster. The two innate 
immune gene clusters, IL-6 signaling/complement and IRG 
cluster, were strongly expressed in ferret lung tissue from 2— 
14 DPI with SARS-CoV relative to the mock primary infection 
group. Conversely, an adaptive immune response gene cluster 
enriched with genes involved in antigen processing and 
presentation was highly expressed in die lungs of SARS-CoV 
infected ferrets after 14 DPI. SARS-CoV titres increased in the 
lungs prior to the initial peak in neutralizing antibody titers at 
14 DPI and then bordered on the limit of detection (Figure 2). 
Upon reinfection with SARS-CoV, the ferrets at 29 DPI did not 
increase gene expression despite an anamnestic response in 
neutralizing antibodies at 3 DPR (Figure 1). We analyzed the 
three clusters of functionally related innate immune genes in 
more detail below. 


Differential Evolution of Gene Expression in Ferret Lungs 
during SARS-CoV Infection and Reinfection 

Our above analysis of the gene expression profiles during 
SARS-CoV infection and reinfection showed a decrease in IRG 
expression during SARS-CoV-reinfection compared to primary 
infection. We next went on to further analyse the expression 
profiles of SARS infection compared to reinfection using the 
DAVID bioinformatic tool [26] was used to perform functional 
classification in differentially expressed genes. During SARS- 
CoV primary ferret infection, a robust increase in the number of 
up- and down-regulated genes on day 2 after infection occurred. 
By day 5 the differential gene expression had decreased 
markedly although it is still moderate levels and on day 7, the 
gene expression was found to be close to baseline. During 
reinfection, a moderate increase in gene expression occurred 
which was maintained for a longer period of time (Table 1). 
Next, we analyzed the pathway activation profiles on 2 DPI 
compared to 2 DPR, which corresponded to peak immune 
activity during both infections (Table 2). On 2 DPI, the gene 
expression patterns showed strong activation of inflammation- 
related genes, activation of MAPK, JAK-STAT and Wnt 
signaling pathways (Table 2). Together, these pathways 
participate in leukocyte activation and migration to infection 
sites. Interestingly, among the upregulated genes on 2DPR, 
lysosomal degradation was the only among the KEGG pathways 
to be significantly enriched (Fisher’s exact test’s p = 1.6E—8); this 
scenario is suggestive of the presence of relevant phagocytosis- 
mediated immunity but without triggering a significant inflam¬ 
matory response. Taken together, these results suggest that there 
are differential gene expression profiles corresponding to SARS 
infection compared to reinfection that may reflect protective host 
immune responses. 
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Figure 2. Viral burden in lung tissue. Viral burdens in sections of lung were determined by the TCID 50 method, as outlined in the materials and 
methods section. LoglO virus titres are shown on the vertical axis vs. study day on the horizontal axis. Note that ferrets were infected on study day 1, 
and reinfected ferrets were also inoculated on study day 30. Mock infected animals received an intranasal instillation of serum-free media on study 
day 1 and had undetectable virus. Values shown represent group means of 3-4 ferrets per group, and error bars show standard deviation. 
doi:10.1371/journal.pone.0045842.g002 


Primary but not Secondary SARS-CoV-infection Increases 
Complement and IL-6 signaling Genes 

The complement cascade plays an integral part in innate 
immunity by labeling pathogens for destruction, inducing leuko¬ 


cyte migration and lysing bacterial [27]. We found that a cluster 
enriched in complement system genes, including C1NH, C1QA/ 
B/C, C1R/S, C3, C4A, C8B, CFB, CR1, FCN1 and MASP1 (see 
Table 3 for full gene names), was significantly upregulated during 
the first 2 weeks post infection (Figure 5A). Quantitative RT-PCR 
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Figure 3. Lung histopathology to SARS-CoV challenge following reinfection. Histological lung sections (5 pm) were obtained from multiple 
lung lobes at 7 d postchallenge and stained by hematoxylin and eosin. Representative micrographs from uninfected (A), SARS infected alone (B), or 
SARS infected and re-infected (C and D) are shown. Primary SARS-CoV infection produced inflammation and the appearance of lung immune cells 
primarily surrounding small-to-medium bronchial airways at 7 days following challenge (B). Ferrets that received infection-reinfection were largely 
protected from lung histopathology (C). 
doi:10.1371 /journal, pone.0045842.g003 
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Figure 4. Microarray analysis of gene expression in lung tissue from ferrets infected and reinfected with SARS-CoV. EDGE analysis 
across all time points identified 3454 genes as significantly differently expressed (>2-fold change in at least one time point, p£0.05, and q<0.1) as 
described in the Methods. Genes were then one-way hierarchically clustered by gene using Pearson correlation and average distance metrics (red, 
upregulated; blue, downregulated). The most significant canonical signaling pathways according to IPA for the resulting clusters are noted. Full gene 
lists are publically available on GEO (see Methods). DPI, days post-infection. DPR, days post-reinfection. 
doi:10.1371/journal.pone.0045842.g004 


(QRT-PCR) validation of C1NH and FCN1 is shown in Figure 
SI. Interestingly, the complement gene cluster was not as 
significantly upregulated the first 2 weeks post-reinfection and 
was subsequently markedly downregulated following 7 DPR. The 
genes of all three complement cascades were integrated at 2 DPI 
determined by Ingenuity Pathway Analysis (IPA), including the 
regulators C1NH and CR1 (Figure 5B). In summary, comple¬ 
ment activation was correlated to the transition from innate to 
specific immune responses and clearance of SARS-CoV during 
primary infection and was not strongly regulated during 
reinfection. 

The pro-inflammatory cytokine, IL-6 is known to function in 
various systems such as B cell development, T cell activation and 
macrophage proliferation [28]. Furthermore, it is a significant 
regulator of fever during infection. We found that the expression of 
IL-6 and seven IL-6 signaling-associated genes, including GRB2, 
JUN, NFKBIA, RRAS2, STAT3, and two MAP kinases (Table 
SI), were increased in ferret lungs as early as 2 DPI and until 
7 DPI (Figure 4C). These genes were not increased following 
reinfection or were decreased. IL-6 and STAT3 gene expression 
was confirmed by QRT-PCR (Figure SI). The context of these 
upregulated molecules at their 5 DPI peak in the IL-6, IL-1, and 
TNF-oc signaling pathways was analyzed by IPA and drawn into 
the respective pathways (Figure 5D). 

Differential IFN Responses in SARS-CoV Infected versus 
Reinfected Ferrets 

In our previous analysis of host immune correlates with 
pathogenic potential in SARS-infected individuals, we reported 
that hyper innate IFN and IFN response gene (IRG) activity could 
be identified in acutely-infected SARS-CoV patients, the persis¬ 
tence of which correlated with ineffective development of adaptive 
immunity and severe clinical course (Table SI) [5,9]. Here we 
found fifty IRGs were significantly upregulated in ferret lung tissue 
during the first 2-7 DPI with SARS-CoV relative to mock- 
infected ferrets, including CD274, IFI30, IFI44, IFI44L, IFI6, 
IRF1, IRF2, ISG15, MNDA, OAS2, and PSMB8 genes (Table 3 


and Figure 6A). QRT-PCR validation was performed on CD274, 
IRF1, IRF2, and MNDA (Figure SI). 

ST ATI and JAK1, key upstream mediators of integrated IRG 
gene expression upon phosphorylation, were significantly upregu¬ 
lated in ferret lungs only during acute SARS-CoV infection. 
STAT1 gene expression was confirmed by QRT-PCR (Figure 
SI). 

Two IFN-induced chemokines, CCL2 and CXCL10, were 
significantly upregulated in the lung at 5-7 DPI with SARS-CoV 
but not after reinfection (Table 3 and Figure 6A). ORT-PCR 
validation of CXCL10 gene expression is shown (Figure SI). 

Interestingly, an additional EDGE analysis of the 3454 genes 
identified above as significantly differentially expressed across all 
time points performed using only the 28 DPI through 28 DPR 
time points identified only 29 significantly increased genes relative 
to 28 DPI, one of which was an IRG, namely ISG15 (Table SI). 

The integrated expression of innate immune genes, IRGs, and 
IFN-mediated chemokines in the lungs of viremic SARS-CoV 
infected ferrets therefore reflected the antiviral responses that 
correlate with only the acute phase of SARS. 

Select IRG Expression in Adjuvant-vaccinated SARS-CoV- 
challenged Ferrets 

The ferret also shows promise in the evaluation of candidate 
SARS vaccines [5,9,29]. Past attempts at developing a CoV 
vaccine have been unsuccessful as many of the candidate vaccines 
caused disease exacerbation, due to cellular or humoral immune 
enhancement via antibodies to the outer envelope protein [9,10]. 
Above, we have shown an immune profile correlating with 
primary SARS-CoV infection in ferrets. Furthermore, SARS-CoV 
reinfection 30 days following primary infection showed a 
significant decrease in innate and adaptive immune gene activity 
associated with immune protection. Therefore as a comparison, 
we next profiled the immune response of SARS-CoV vaccinated 
ferrets challenged with SARS-CoV to determine if a protection 
molecular signature could be induced following vaccination. 

To investigate the a vaccination molecular signature, ferrets 
were mock-vaccinated as control or vaccinated using an inacti- 


Table 1 . Number of regulated genes in different functional categories following primary infection and reinfection. 




Days post primary infection 

Days post reinfection 




2 

5 

7 

2 

5 

7 

Total 

T 

3025 

1091 

231 

1151 

819 

1773 


1 

1348 

449 

107 

228 

534 

473 

Cellular process 

T 

1779 

668 

144 

646 

444 

1028 


1 

683 

245 

57 

117 

281 

255 

Metabolic process 

T 

1296 

488 

95 

469 

307 

754 


1 

489 

175 

47 

83 

204 

181 

Response to stimulus 

T 

483 

179 

66 

173 

128 

273 


1 

221 

85 

17 

44 

93 

84 


Number of regulated genes in different functional categories with at least 1.5-fold change and a significant t-test of p<0.05 (f upregulated, 1 downregulated). 
doi:10.1371/journal.pone.0045842.t001 
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Table 2. Intersect analysis of upregulated genes in the lungs of ferrets after infection or reinfection with SARS-CoV. 




2 DPI only 

2 DPI & 2 DPR 

2 DPR only 

Total upregulated genes 

2626 

399 

752 

Gene Ontology 

Signal transduction 

321 

41 

85 

Immune system process 

121 

11 

26 

Ubiquitin-dependent protein catabolic process 

64 

6 

16 

G-protein coupled receptor protein signalling pathway 

63 

4 

17 

Protein kinase cascade 

53 

5 

11 

Inflammatory response 

45 

6 

10 

KEGG Pathways 

MAPK signaling pathway 

44 

5 

7 

Wnt signaling pathway 

26 

4 

3 

Cytokine-cytokine receptor interaction 

24 

0 

4 

Jak-STAT signaling pathway 

19 

3 

4 

Cell adhesion molecules (CAMs) 

18 

3 

4 

Lysosome 

14 

9 

18 

Toll-like receptor signaling pathway 

13 

2 

1 


"2DPI only": specific of 2 days after infection, "2DPR only": specific of 2 days after reinfection, and "2 DPI & 2DPR": genes upregulated in both experimental groups. CoV 
(^1.5-fold change and Student's t-test p<0.05). 
doi:10.1371/journal, pone.0045842.t002 


vated whole virus SARS-CoV vaccine candidate formulated with 
or without Al(OH) 3 , SARS FCP-A and SARS FCP respectively 
[30]. EDGE microarray analysis was then preformed on lung 
RNA collected at 2, 7, and 28 days post SARS-CoV challenge 
(DPC). Of the 3717 significantly altered genes resulting from this 
combined time course analysis, 42 were identified as IRGs by IPA 
(Figure 6A), many of which had also been identified in the 
infection-reinfection experiment (Table 3). When we applied a 
two-fold mean difference biological filter between the three groups 
at each time point, however, it was found that 8 of the 42 IRGs 
were downregulated at least two-fold at 2 or 7 DPC in ferrets 
previously immunized with SARS FCP-A relative to the mock- 
vaccinated group (Figure 6B). These 8 IRGs, included: ISG15 
(QRT-PCR shown in Figure SI), CCL2, IFI44, and PSMB8 
discussed above; GPX1 (an antioxidant enzyme), and PSME1 and 
2 (roles in MFIC class I antigen processing, see Table 3 for full 
gene names); and myxovirus resistance 1 (MX1, a potent antiviral). 
Interestingly, the temporal expression patterns of 4 IRGs (ISG15, 
IFI44, PSME2 and CCL2) were similar amongst the SARS-CoV 
infection-reinfection experiments (50 IRGs) and this list of 8 IRGs 
from the adjuvanted vaccine-challenge experiments (Figure 6C). 
These data show a downregulation of early SARS-CoV challenge 
regulated ISG genes in the context of a SARS vaccine. 

Discussion 

Currently there exists a paucity of information surrounding the 
molecular events associated with protection from SARS-CoV 
infection. Furthermore, due to the evolutionary and promiscuous 
nature of the SARS-CoV and other coronaviruses there is a 
current need to develop vaccination and therapeutic strategies for 
humans and mammals. Our time course analysis of differential 
gene expression in the lungs of SARS-CoV infected and SARS- 
CoV reinfected ferrets identified three key clusters of functionally 
related genes. Early (2-14 DPI) expression of IL-6 signaling/ 
complement and IFN response genes followed immediately 


(14 DPI) by antigen processing and presentation gene expression 
correlated with peak SARS-CoV titres in the lungs and peak 
neutralizing antibody titers, respectively. Reinfecting the ferrets 
with SARS-CoV did not reinitiate the same organized expression 
of antiviral innate immune response genes; however adaptive 
immunity in the form of SARS-specific antibody production 
quickly appeared at 3 DPR. These correlations suggested that a 
period of IFN-driven innate antiviral responses mediates acute 
SARS-CoV infection, after which specific adaptive immune 
responses confer protection. These findings are specifically 
imperative to the understanding of SARS-CoV and SARS-like- 
CoV infections and provide a comprehensive baseline for the 
development of CoV antiviral therapeutics and vaccines tailored 
to account for the specific IFN responses and subsequent antibody 
production. 

A protective role for complement in SARS-CoV pathogenesis 
has been proposed [31]; however, the extent to which complement 
is involved in SARS-CoV infection is unknown. In this study we 
found that certain complement system genes were significantly 
upregulated during the first 2 weeks post infection. CFB is key to 
the alternative complement cascade and is cleaved into Ba and Bb 
(UNIPROT). The Bb serine protease combines with C3b to 
generate the C3 convertase. Clq is the first target recognition 
protein of the classical complement cascade and is an important 
link between innate and adaptive immunity [32]. Cl complex 
activity is regulated by C1NH which forms a complex with Clr, 
Cls, and MASP proteases. We have previously shown that C1NH 
is upregulated in SARS-CoV infected patients prior to either 
recovery or progression to severe disease [33], C4 is cleaved by 
activated Cl to produce C4b, an essential component of the C3 
and C5 convertases of the classical pathway (UNIPROT). FCN1 
initiates the lectin complement cascade via MASP zymogens, e.g. 
MASP1, which in turn cleave C4 [34], Interestingly, the 
complement signaling regulators, C1NH and CR1, were also 
upregulated at 2 DPI. 


PLOS ONE | www.plosone.org 


7 


September 2012 | Volume 7 | Issue 9 | e45842 
















IFN Responses of SARS CoV Infection in Ferrets 


(A) Interferon Response Genes 



-4 0 4 

Fold Change (Log 2 ) 

(B) Interferon Signaling Pathway Analysis 
3 DPI 



o 

Complex 


Cytotdne/Growth Factor 


Enzyme 


Kinase 


Peptidase 

H 

Phosphatase 

cc 

Transcription Regulator 

V 

Transmembrane Receptor 

o 

Transporter 

o 

UnKnown 


■Relation 


-4 0 4 

Fold Change (Log 2 ) 


$ S xtu* OAS2 mnda ISG15 $ 


PLOS ONE | www.plosone.org 


8 


September 2012 | Volume 7 | Issue 9 | e45842 





































IFN Responses of SARS CoV Infection in Ferrets 


Figure 5. Complement and IL-6 signaling in SARS-CoV infected-reinfected ferret lungs. (A) Complement and IL-6 signaling (C) genes 
selected by pathway analysis are shown in a one-way hierarchical cluster (red = upregulated, blue = downregulated). (B) IPA canonical complement 
system pathway analysis at 2 DPI. (D) IPA canonical IL-6 signaling pathway analysis at 5 DPI. All genes are significantly differently expressed (EDGE 
analysis: >2-fold change in at least one time point, p£0.05, and qsO.1). 
doi:10.1371 /journal.pone.0045842.g005 


IL-6 has been shown to be induced in human bronchial 
epithelial cells following SARS-CoV infection in vitro and in vivo 
as well as associated with hyper-immune activation during 
SARS-CoV infection [14,35,36]. We found that the expression 
of IL-6 and IL-6 signaling-associated genes were increased in 
ferret lungs for a longer period of time than complement genes 
(2 DPI to 7 DPI). IL-6-regulated genes are mainly induced by 
IL-6/IL6 receptor signaling via JAK2:STAT3 and the RAS/ 
MAP/NF-IL6 pathway [37]. IL-6 gene expression is also 
partially dependent on NF-kB which is activated by NFKBIA 
(I-kBo() degradation [38]. Since IL-6 signaling can drive further 
IL-6 expression via the RAS/MAP kinase cascade and NF-IL6 
activation [37] our results suggest that potentially self-sustaining 
acute phase responses occur during early SARS-CoV infection 
and not after reinfection. 

Critical to IFN-mediated antiviral activity is JAK/STAT 
signaling which prompts widespread transcriptional activation of 
IRGs [39], Resembling the results of our PBMC study of acute 
SARS-CoV infection [33] and the previous results of a SARS- 
CoV mouse model [14,40], a prominent number of IRGs were 
upregulated in ferret lung tissue during the first 2-7 DPI with 
SARS-CoV infection relative to mock-infected ferrets. CD274 is 
induced by IFNs and functions in T cell costimulation during viral 
infection [41], IFI30 and PSMB8 are induced by IFN-y and have 
roles in antigen presentation (UNIPROT). IFI44, IFI44L, IFI6, 
MNDA, and OAS2 are classical antiviral IRGs induced mainly by 
type I IFNs (UNIPROT). IRF1 is rapidly induced by IFNs and 
binds to the upstream regulatory region of type I IFN and IFN- 
inducible MHC class I genes [42]. IRF2, conversely, binds to the 
interferon consensus sequence and represses IRF-mediated genes 

[43]. 

While we show that IRF/IRG gene expression is integrated 
during acute infection in our non-fatal ferret SARS-CoV 
infection model, we previously identified dysregulated IFN 
activity in SARS-CoV infected patients with severe clinical 
course [33], suggesting that differently balanced IRG gene 
signatures affect antiviral versus proinflammatory host immune 
response development [44], STAT1 phosphorylation and IRG 
transcription invokes a cellular antiviral state where STAT1 has 
been previously reported up-regulated during SARS-CoV 
infection [13,45]. STAT1 and JAK1 were indeed significantly 
upregulated in ferret lungs only during acute SARS-CoV 
infection. STAT1-deficient mice suffer worse disease upon 
SARS-CoV infection with greater systemic infection and viral 
burden in the lungs [13,46]. 

Two IFN-induced chemokines, CCL2 and CXCL10, were 
significantly upregulated in the lung at 5-7 DPI with SARS-CoV 
infection but not after reinfection. CCL2 protein expression has 
been associated with early inflammatory host responses during 
murine SARS-CoV infection [47]. We previously demonstrated 
that high levels of CXCL10 protein were present in plasma from 
SARS-CoV infected patients during early illness and quickly 
returned to baseline levels in those patients that quickly 
recovered [33,48]. In severe cases of SARS-CoV infection, 
however, CXCL10 levels remained significantly elevated for the 
duration of the patient’s infection. Interestingly, ISG15 was the 
only known IRG that was significantly increased following 
reinfection. ISG15 is an antiviral ubiquitin-like modifier that 
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conjugates with JAK1 and STAT1 after type I IFN stimulation 
[49], Also, antiviral ISG15 derivatives have been shown to be 
preferred substrates for the deubiquitinating activity of the 
SARS-CoV papain-like protease [50]. ISG15 is expressed by a 
wide variety of lymphoid cells and tissues. ISG15 and its targets 
may therefore represent unique innate immunity correlates that 
are not influenced by tissue viral loads, but that may participate 
in the IFN-mediated transition from innate to adaptive 
immunity. 

Collectively, the marked expression of IFN-mediated responses 
in die lungs of viremic SARS-CoV infected ferrets and not 
reinfected ferrets further asserts that proinflammatory IFN 
responses complement the acute phase of SARS and that 
resolution of IRG activity is associated with priming of an 
anamnestic response diat will neutralize SARS-CoV reinfection 
without reinitiating acute inflammation. Although, it is important 
to mention the lack of upregulation following second infection may 
be due to the low virus replication levels after secondary infection. 
The presence of a low immune stimulation in the reinfection 
group is probably associated with the marginal viral replication. 
Also, the analysis of the gene expression profiles is suggestive of a 
scenario with very limited immune stimulation and without 
prominent participation of suppressor genes. 

The general lack of IFN-mediated immune responses in ferrets 
reinfected with SARS-CoV suggested that a SARS vaccine- 
challenge experiment would reveal additional immune correlates 
of immunity. Our microarray analysis on SARS-CoV challenged 
ferret lungs with or without previous immunization using an 
Al(OH) 3 -adjuvanted, inactivated whole virus SARS-CoV vaccine 
candidate revealed die downregulation of 8 IRGs. Moreover, four 
IRGs (ISG15, IFI44, PSME2 and CCL2) were expressed in 
common between the SARS-CoV infection-reinfection experi¬ 
ments and these 8 IRGs. The select abrogated expression of IRGs 
in SARS FCP-A vaccinated, SARS-CoV challenged ferrets refiect 
the infection-reinfection results demonstrating lack of IRG 
induction upon reinfection. The expression of these genes early 
in the natural course of infection suggested IRGs to be indicative 
of an early stage in SARS-CoV disease, but neither reinfection nor 
SARS vaccination-challenge provided evidence for IRGs as 
biomarkers of immune responses that are capable of suppressing 
SARS-CoV challenge. Importantly, the focus of this study was on 
the gene expression and protein analysis would be a valuable 
addition to future experiments. Furthermore, when we measured 
SARS-CoV viral RNA levels in SARS-CoV challenged mock-, 
SARS FCP-, and SARS FCP-A-vaccinated ferrets lungs by PCR, 
all nrock-vaccinated ferrets were PCR-positive for SARS-CoV at 
7 DPI, one SARS FCP-vaccinated ferret was positive at each of 2 
and 7 DPI, and no SARS FCP-A-vaccinated ferrets were PCR- 
positive for SARS-CoV at any time. 

Concerns regarding immune enhancement arose when immu¬ 
nization of ferrets with recombinant modified vaccinia Ankara 
(rMVA) expressing SARS-CoV spike induced strong inflammation 
and hepatitis [29], later attributed to rMVA expressing SARS- 
CoV antigens [9]. More recently, formalin-inactivated whole-virus 
[24] and adenovirus-based [22,23] SARS-CoV vaccine candidates 
have shown promise in reducing pneumonia during SARS-CoV 
challenge in ferrets. Interestingly, immune enhancement was not 
seen in mice vaccinated with Al(OH)3-adjuvanted or non- 
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Table 3. Expression of selected significantly changed genes over time in SARS-CoV infected-reinfected ferret lungs. 
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Figure 6. IFN responses in SARS-CoV infected-reinfected ferret lungs. (A) Fifty IRGs selected by pathway analysis are shown in a one-way 
hierarchical cluster. (B) IPA canonical IFN-signaling pathway analysis at 3 DPI and 3 DPR. All genes are significantly differently expressed (EDGE 
analysis: >2-fold change in at least one time point, p<0.05, and q<0.1). 
doi:10.1371/journal.pone.0045842.g006 


adjuvanted inactivated SARS-CoV whole virus candidate vaccines 
[30]. Furthermore, vaccine strategies for coronaviruses in an aged 
mouse model have also been investigated and found that wild type 
VRP (Venezuelan equine encephalitis virus replicon particle) 
vaccines protected animals from challenge [15]. Here, IFN- 
mediated immune enhancement was not noted in our ferret 
SARS-CoV infection-reinfection model and we have previously 
reported that increased neutralizing antibody titers after SARS- 
CoV reinfection is a correlate of immune protection [51]. One 
caveat, due to the large nature of our study, was the absence of an 
Alum control only group and it would be important in the future 
to investigate the genes regulated by Alum treatment alone. 

Taken together, this study has given insight into the molecular 
events following SARS-CoV infection as well as identified an 
imperative signature for immunological protection against SARS- 
CoV in both a reinfection and vaccination model. Specifically, the 
ferret model of SARS-CoV infection showed a two step temporal 
paradigm of host immune responses in the lung. SARS-CoV 
infection in the ferrets was characterized by an early innate 
immune response (robust IFN and chemokine gene expression) 
followed by adaptive immunity likely eliciting localized anti- 
SARS-CoV antibodies and protective immune responses. Like¬ 
wise, the lack of repeat robust innate immune responses in ferrets 
when reinfected with SARS-CoV suggested that the post¬ 
reinfection gene signatures may reveal immune correlates related 
to successful immunization that may be applied to future use of the 
ferret model of SARS infection in vaccination strategies. Thus, the 
identified IRGs have the potential to be inverse correlates of 
immunity against acute SARS-CoV infection versus reinfection, or 
SARS-CoV vaccination versus challenge. This dataset serves as a 
resource in human SARS-CoV vaccine candidate trials for 
modeling vaccine efficacy in populations at risk for severe 
outcome. 


Methods 

Ethics Statement 

All work with animals was conducted in strict accordance with 
the Institutional Animal Care and Use Committee (IACUC) 
guidelines from Southern Research Institute (SRI) or Lovelace 
Respiratory Research Institute (LLRI). For SRI, ferret studies 
were approved by the Southern Research Institute’s Institutional 
Animal Care and Use Committee. Southern Research Institute 
has Veterinary Medicine tasked to monitor and support all animal 
experiments. Research was conducted in compliance with the 
Animal Welfare Act and other federal statutes and regulations 
relating to animals and experiments involving animals and adheres 
to principles stated in the Guide for the Care and Use of 
Laboratory Animals, National Research Council, 1996. The 
facility where this research will be conducted is fully accredited by 
the Association for Assessment and Accreditation of Laboratory 
Animal Care International. For LLRI, The animal work was 
reviewed and approved by the LRRI Institutional Animal Care 
and Use Committee under protocol number 06-017. All work 
involving infectious agents are reviewed by the LRRI Infectious 
Agents Committee and approved. All infections and sample 
collections were performed under 5% isofluorane anesthesia and 
all efforts were made to minimize suffering. 

Cell Culture, Viruses and Vaccine 

SARS-CoV TOR2 was amplified in Vero E6 cells (ATCC, 
VERO Cl008, CRL 1586) at the Southern Research Institute 
BSL-3/ABSL-3 facility as described. The SARS-CoV strain used 
for vaccine preparation was the CDC strain, AY714217 [30]. 

Study Design 

Influenza-free, castrated and descented male Fitch ferrets 
(.Mustela putorius furo), 36 to 45 weeks of age, were purchased from 
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Triple F Farms (Sayre, PA). Infection-reinfection: 72 ferrets were 
randomly assigned to groups: ( 1 ) mock primary infection (n = 18), ( 2 ) 
primary infection (n =18), (3) mock infection-reinfection (n = 18), 
and (4) infection-reinfection (primary and secondary inoculation) 
(n = 18). On day 0, the group two and four ferrets were anesthetized 
and received 500 pL of serum-free medium containing SARS-CoV 
per nare at a total dose of 10 3 TCID 50 . On 29 DPI, i.e. 0 DPR, 
group four ferrets were inoculated again with 10 3 TCID 50 SARS- 
CoV. The mock groups received 1 ml of serum-free medium. Three 
group one and two ferrets were sacrificed at 2, 3, 5, 7, 14, and 
28 DPI. Three group three and four ferrets were sacrificed at 2, 3, 5, 
7, 14, and 28 DPR. Vaccine—challenge: 57 ferrets were assigned to 
groups: (1) mock vaccine with SARS-CoV challenge (n= 18), (2) 
SARS-CoV vaccine with no adjuvant and SARS-CoV challenge 
(n = 18), (3) SARS vaccine with adjuvant and SARS-CoV challenge 
(n = 18), and (4) mock vaccine-challenge (n = 3). Vaccinated ferrets 
received a double-inactivated whole virus SARS-CoV vaccine 
candidate (Baxter Innovations GmbH, Vienna, Austria, Bulk Drug 
Substance: SA/0/05/01/R/PBDS) formulated with (SARS FCP- 
A, 5pg/ferret/dose) or without (SARS FCP, 10 pg/ferret/dose) 
Al(OH) ;i as adjuvant [30]. Mock-vaccinated ferrets received 
TBS+0.05% poloxamer (Baxter, Product Specification ^SP0050). 
Inoculation with 10 3 TCID 50 SARS-CoV occurred 42 days post 
primary vaccination with a boost delivered two weeks prior to 
infection. Six ferrets from groups one, two, and three were sacrificed 
at 2, 7, and 28 DPC. Mock vaccine-challenge ferrets were 
sacrificed at 2 DPC. All procedures were in accordance with the 
NRC Guide for the Care and Use of Laboratory Animals, the 
Animal Welfare Act, and the CDC/NIH Biosafety in Microbio¬ 
logical and Biomedical Laboratories. All animal experiments were 
conducted in a registered BSL-3/ABSL-3 facility at the Lovelace 
Respiratory Research Institute. 

Viral Neutralization Assay and Quantification 

Neutralization tests in serum and TCID 50 quantification of 
SARS-CoV in lung tissue suspensions were performed as 
described. Virus titres were expressed as TCID 50 units (U)/g. 
The lower limit of detection was 10 2 TCID 50 U/g. RT-PCR 
detection (within 30 cycles) of SARS-CoV was performed on 
isolated lung RNA as previously described [52] using ORFlb F:5'- 
CAGAACGCTGTAGCTTCAAAAATCT-3', R:5 '-TCA- 

GAACCCTGTGATGAATCAACAG-3' primers. 

Microarray Analysis 

Total RNA from lung was prepared for microarray analysis as 
described, as well as microarray dataset preprocessing and 
normalization against the appropriate mock dataset [53]. RNA 
quality from one infection-reinfection group four ferret’s 3 DPR 
lung was insufficient for arraying. Affymetrix Genechip Canine 
Genome 2.0 arrays (Affymetrix, Santa Clara, CA) were used to 
assess ferret-like gene expression as previously established [53,54], 
Infection-reinfection: EDGE [25] was used to identify significantly 
differentially expressed genes over all time points for groups two 
and four (2 DPI -28 DPR) using the 2 DPI mocks as baseline. 
Vaccine-challenge: EDGE was used to identify significantly 
differentially expressed genes in each of groups one, two, and 
three using the 2 DPC group 4 datasets as baseline, the results of 
which were combined. Hierarchical clustering, functional pathway 
analysis (IPA, Ingenuity Systems, Redwood City, CA), and gene 
annotation (includes the Universal Protein Resource (UNIPROT, 
URL: www.expasy.org/uniprot) was performed as previously 
described [53]. Datasets are available at the NCBI Gene 
Expression Omnibus, ^GSE11704. 


QRT-PCR 

QRT-PCR was performed in triplicate as described [53]. Ferret 
primer sequences: C1NH, F:5'-GCCTCTCAGAGCCTGTAT- 
GG-3', R:5'-CTTCCACTTGGCACTCAGGT-3', CD274, 
F:5 '-GGCAATGTGACAATGGAATG-3', R:5 '-CTCTGGC- 
TGTAGCTGCTGTG-3', CXCL10, F:5'-CTTTGAACCAAA- 
GTGCTGTTCTTATC-3', R:5'-AGCGTGTAGTTCTAGA- 
GAGAGGTACTC-3', FCN1, F:5'-CACCAAGGACCAGGA- 
CAATGA-3', R:5'-CACCAGGCCCCCTGGTA-3', IL6, F:5'- 
AGTGGCTGAAAC AC GTAAC AATTC-3', R: 5'-ATGGCC- 
CTCAGGCTGAACT-3', IRF1, F:5'-CTACCTGCAGCCTC- 
CACTTC-3', R:5'-GAGACCAACACGGTCAGGTT-3', IRF2, 
F: 5' - AGGT G AC C AC C G AG AGT G AC - 3', R:5'-CCCCATG- 
TTGCTGAGGTACT-3', ISG15, F:5'-AGC AGC AGATAGC- 
CCTGAAA-3', R:5 '-C AGTTCTTC ACC ACC AGCAG-3', 

MNDA, F:5'-CCTGCGTGGACAAGCTAATA-3', R:5'- 
GCTGCTTCTTCTTGGGAGTT-3', STAT1, F:5'-AGCCTT- 
GCATGCCAACTCA-3', R:5'-ACAGTCCAGCTTCACCGT- 
GAA-3', and STAT3, F:5'-CAACCCCAAGAACGTGAACT- 
3', R:5'-AGCCCACGTAATCTGACACC-3'. 

Histopathology 

On Day 7 pi, animals were euthanized and the upper and lower left 
lobes of the lung tissues were harvested and formalin-fixed by perfusion 
followed by paraffin embedding and sectioning. Tissue slides were then 
stained with hematoxylin and eosin for histopathology assessment. 

Statistical Analysis 

The microarray analysis used EDGE time course modeling [54] 
with p<0.05 and q£0.1 for significance and false discovery rate, 
respectively, and &2-fold change in at least one time point. 

Supporting Information 

Figure SI QRT-PCR Validation of microarray analysis of 
ferret lungs infected and reinfected with SARS-CoV or 
vaccinated and SARS-CoV infected. QRT-PCR analysis of gene 
expression in lung tissue from ferrets infected and reinfected with 
SARS-CoV (AJ). QRT-PCR analysis of ISG15 expression in lung 
tissue front SARS vaccinated-challenged ferrets. Values shown 
represent the groups’ mean mRNA levels for the indicated time 
points relative to the mean mRNA levels for the appropriate mock 
animals at 2 DPI (A-J) or 2 DPC (K). All analyses were performed in 
triplicate. Group means represent data from 3 ferrets. DPI, days post¬ 
infection. DPR, days post-reinfection. DPC, days post-challenge. 
Error bars represent standard error. * p = 0.002 (Student’s t-test) 
difference from the mock, ** non-significant difference from the mock. 
(TIF) 

Table SI Genes significantly changed over time in 
SARS-CoV infected ferret lungs following reinfection. 

(DOC) 
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